PQCD analysis of exclusive e+e -^ D'^.D, , cross section at ^i = 10.6 GeV 
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We analyze the exclusive pseudoscalar Df,D7, pair production in e^e~ annihilations at -^s — 10.6 
GeV using a non-factorized PQCD with the light-front wave function that goes beyond the peaking 
approximation. We compare our non-factorized analysis with the usual factorized analysis based 
on the peaking approximation in the calculation of the cross section for the heavy meson pair 
production. We also discuss the higher helicity contribution to the cross section. Our analysis 
provides a constraint on the size of quark transverse momentum inside the D meson from the recent 
Belle data, (TExp.(e+e~ -^ D+D~) < 0.04. 



I. INTRODUCTION 

Recently, many theoretical works have been devoted 
to explain the large discrepancy between theoretical and 
experimental results for the charmonium production in 
e+e" annihilations. For instance, the data 0, 0, y] for 
charmonium production cross sections in e~^e~ annihi- 
lations at the S-factory energy i/s — 10.6 GeV differ 
a lot from the theoretical predictions for both exclu- 
sive II, IM B 13 and inclusive Q processes although 
higher order corrections may reduce the differences [lOJ. 
A couple of years ago, the Belle Collabotation j3j re- 
ported the first measurement of the e+e^ — > D*'^D*^ 
and e~^e~ —^ D~^ D*~ cross sections and polarizations at 
^/s ~ 10.6 GeV. They also set an upper limit on the 
cross section for e+e^ —» D^D^. Interestingly, while 
the theoretical predictions based on the heavy quark ef- 
fective theory 4] and the constituent quark model j3| for 
e+e^ -^ D*^D*^ and e+e^ — > D^D*^ cross sections 
are similar to the measured data i^j, the predictions for 
e+e" -^ D^D^ cross section are either quite smaller Q 
or somewhat larger [3 than the data Q- 

The above exclusive/inclusive meson pair productions 
provide a unique opportunity to investigate asymptotic 
behaviors of various meson form factors in the frame- 
work of perturbative quantum chromodynamics(PQCD). 
The heavy meson pair production is of special interest 
since gluons carrying large momentum transfers can be 
rather easily accessible in the kinematic region above the 
threshold. Also, the wavefunctions of heavy systems may 
be well constrained due to the heaviness of constituents. 
Thus, it has been pointed out that exclusive pair pro- 
duction of heavy mesons can be reliably predicted within 
PQCD [ij. 

If the factorization theorem in PQCD is applicable to 
exclusive processes, then the invariant amplitude for ex- 
clusive processes factorizes into the convolution of the 
valence quark distribution amplitude(DA) (j)(x,q^) with 
the hard scattering amplitude Tr [131 ■ To implement the 
factorization theorem at high momentum transfer, the 
hadronic wave function plays an important role linking 
between long distance nonperturbative QCD and short 
distance PQCD. A particularly convenient and intuitive 



framework in applying PQCD to exclusive processes is 
based upon the light-front (LP) Fock-state decomposi- 
tion of hadronic state. In the LP framework, the va- 
lence quark DA is computed from the valence LP wave 
function \['„(xi,k^i) of the hadron at equal LP time 
T = t + z/c which is the probability amplitude to find 
n constituents(quarks,antiquarks, and gluons) with LP 
momenta ki = (a;i,k_Li) in a hadron. Here, Xi and kj_i 
are the LP momentum fraction and the transverse mo- 
menta of the ith constituent in the n-particle Pock-state, 
respectively. 

To lowest order in perturbation theory of the meson 
form factor calculation at large momentum transfers, the 
hard scattering amplitude Th is dominated by one-gluon 
exchange diagrams. Por the factorization theorem to be 
applicable in the heavy meson pair production analysis, 
the only consistent form of the quark DA would be the 5 
function, i.e. <j){x, q^) ~ 5{x — ttiq/M) where tuq and M 
are the heavy quark mass and the meson mass, respec- 
tively 13]. In this so called "peaking approximation", 
the momentum fraction carried out by ith constituent 
is equal to the ratio of the constituent mass to meson 
mass, Xi = vfiijM. This relation implies M — mi 4- 7712, 
i.e. zero-binding energy limit. 

However, as pointed out in Ref. [l^l, if the quark DA 
is not an exact 6 function, i.e. k^ in the soft bound 
state LP wave function can play a significant role, the 
factorization theorem is no longer applicable. To go be- 
yond the peaking approximation, the invariant amplitude 
should be expressed in terms of the LP wave function 
'i{xi,'k±i) rather than the quark DA. In Ref. J^], the 
validity issue of peaking approximation for the heavy me- 
son pair production processes was discussed using the LP 
model wave function '^{xi,'k±i) ex exp(— Mq//3^), where 
Mq is the invariant mass of the constituent quark and 
antiquark defined by Mq = J2i(M-\ + ''^D/^i ^.nd (3 is 
the gaussian parameter. The limit /3 — > corresponds 
to the peaking approximation(i.e. zero-binding energy 
limit M = nil + 7712). In the analysis of the heavy-heavy 
system like Bc{bc) meson, it was found that the effect of 
going beyond the peaking approximation (/3 up to 100 
MeV) was not important cornmred to the peaking ap- 
proximation limit(i.e. /3 ^- 0) [I3. However, it is not yet 
clear if the same conclusion would apply to the heavy- 



light system such as D and B mesons. Moreover, the 
initial analysis limited only up to /3 < 100 MeV may not 
be sufficient to draw a definite conclusion on the validity 
of the peaking approximation. 

The main purpose of this work is to extend the previous 
analysis [13 and point out that the recent Belle data Q 
can provide a rather stringent constraint on how broad 
or narrow the D^ meson quark DA is. Clarifying the 
relation between the (3 value and transverse momentum, 
k_L, is a particularly important issue since the quark DA 
is very sensitive to the (3 value and the different shape of 
the quark DA could enhance or reduce the cross section 
for the exlcusive meson pair production in e^e~ annihila- 
tions. Incidentally, Bondar and Chernyak 14] considered 
a rather broad quark DA(i.e. rather significant binding 
energy effect) instead of (5-type quark DA to explain the 
data for the exclusive e+e^ -^ J/t/j + ric process. Ma and 
Si '15] also previously discussed the variation of DA to 
explain the data for the same process. 

In this work, we stress a consistency of our analysis 
in going beyond the peaking approximation. In particu- 
lar, we confirm that the /3 value in our model LF wave 
function is related with the transverse momentum via 
f^QQ ~ V0^'±)oo' ^^ expected, the non-zero /3 value cor- 
responds to the transverse size of the meson and /3 — > 
limit corresponds to the peaking approximation(i.e. zero 
binding energy limit) as discussed in [I^- This implies 
that it may be significant to keep the transverse momen- 
tum kj^ both in the wavcfunction part and the hard scat- 
tering part together before doing any integration in the 
amplitude if f3 is not so close to zero or the binding en- 
ergy effect is not negligible. Thus, we think that the 
factorization of amplitude by integrating out the trans- 
verse momentum seperately in the wavefunction part and 
in the hard scattering part may not provide a consistent 
analysis to take into account the binding energy effect. 
This could distinguish our method from Ref. 14] to take 
into account the binding energy effect. 

We also note that our gaussian parameter /3 is not 
chosen arbitrarily but fixed by the variational principle 
for the well-known linear plus Coulomb interaction mo- 
tivated by QCD llq . which in turn uniquely determine 
the shape of the quark DA in our model calculation. This 
implies that the recent data by the Belle collaboration Q 
provide a useful test on our model calculation. 

The paper is organized as follows. In Sec. II, we 
describe the formulation of our light-front quark model 
(LFQM), which has been quite successful in describing 
the static and non-static properties of the low- lying 
mesons |l6l Il7|. In Sec. Ill, the transverse momentum 
dependent hard scattering amplitude for the meson is 
given within the LF framework. The contribution to the 
meson form factor from higher hclicity components is also 
given in this section. In Sec. IV, the analytic continuation 
from spacelike region to the timelike region is introduced 
to obtain the cross section for the pseudoscalar meson 
pair(Af Af ) production in e~^e~ annihilations. As a valid- 
ity check of our model, we also show that our result for 



the meson form factor obtained in Sec. Ill reduces to the 
peaking approximation in the /3 — > G(i.e. zero-binding) 
limit. In Sec. V, we present the numerical results for 
the e+e^ -^ I?7J,Z3r, cross section and compare with the 
available data. Summary and conclusions follow in Sec. 
VI. In the Appendix, we briefly summarize our proof of 
vanishing contribution from the light-front gauge part in 
M = Mo limit. 



II. MODEL DESCRIPTION 



In our LFQM, the meson wave function is given by 

*5[/^(x,ki,AA) = (^(a;,ki)7^5[(^(x,ki). (1) 

where (/>(x,k_L) is the radial wave function and 
TZ'^'^' {x,'k±) is the spin-orbit wave function obtained by 
the interaction-independent Melosh transformation [13] 
from the ordinary equal-time static spin-orbit wave func- 
tion assigned by the quantum numbers J^^ . The meson 
wave function in Eq. (Q is represented by the Lorentz- 
invariant variables Xi = pf /P~^ , ]s.±i = p±i — XiP± and 
Ai, where P,Pi and A^ are the meson momentum, the 
momentum and the helicity of the constituent quarks, 
respectively. 

The radial wave function (f>(x,'k±) of a ground state 
pseudoscalar meson (J^*^ = ^) is given by 



(a;,k^) = 



7,3/2^3 



1/2 



exp(-fcV2/32), (2) 



where the gaussian parameter P is related with the size 
of the meson. Here, the longitudinal component kz of the 
three momentum is given by k^ — {xi — ^)Mq + (1712 — 
m1)/2MQ with the invariant mass 



M' 



ki 



ki 



Xl 



X2 



(3) 



where Xi = x and X2 = 1 — x. The covariant form of the 
spin-orbit wave function Tl^°^{x, k±) for the pseudoscalar 
meson is given by 






w(pi,A)75w(p2,A) 



V2[M2-(mi-m2)2]i/2' 
and its explicit matrix form is given by 



(4) 






-k^ 
-xim2 — X2mi 



Xim2+X2mi ■ ,^s 

i,R . (5) 



where C = y/2xiX2[M^ - (mi - 7^2)^] and fc^(^) = k^± 

iky. Note that X^aa -^aa aX ~ ^' "^^^ normalization of 
our wave function is given by 

^ fd'k\^>f,{x,-k^,XX)\^ 
XX •' 

= J\x J d'-k^(^^ycl,{x,kj_)\' = 1, (6) 



where the Jacobian of the variable transformation 
{x, k±} -^k = {'k±,k^) is given by 



dk. 



Mo 



dx 4a;ia;2 



(mi — 7712)^ 



M2 



(7) 



The effect of the Jacobi factor has been analyzed in 
Ref. iiaj- 

With this normalization, the root-mean-square(r.m.s.) 



value of the transverse momcntum(i /(k^)QQ) is ob- 
tained via 



(ki>QQ = eg / d^fc|ki||0(x,k^)|^ + (Q^Q). (8) 



Numerically, we confirm that J (k\)QQ = Pqq. The 

numerical values of Pqq are discussed in Sec. V (see 
Table I). 

The quark distribution amplitude(DA) of a meson, 
4>M,x{x, Q), i.e. the probability of finding coUinear quarks 
up to the scale Q in the L^ = O(s-wave) projection of the 
meson wavefunction |l^ is defined by 



4>M,\{x,Q) 



Q 



[d^k^]^{x,-k^,X~X), 



(9) 



where [(fk_i_] = (fkj_y/dk^/dx/Vl6^ for ^ = *55. 



III. HARD SCATTERING AMPLITUDE WITH 
kx-DEPENDENCE 

In this section, we calculate the pseudoscalar me- 
son electromagnetic form factor in the region where the 
PQCD is applicable. Our calculation is carried out using 
the Drell-Yan-West frame ^{q+ = g° + g^ = 0) with 
q^ = Q^ — —(f'- The momentum assignment in the 
q^ — Q frame is given by 

P^(P^|;,o.),P'^(P^^^,q.) 



q = (0,^,q^), 



(10) 



where prime denotes the final state momentum and q = 
P' ~ P and M is the physical meson mass. 

As a starting point, the electromagnetic form factor of 
a pseudoscalar meson is given by a convolution of initial 
and final meson wavefunctions: 



xM'50(x,ki,AA), (11) 



where cfik = dx(Pk±y/dkJdx, k'_L = k^ +2^2^^ and 
is the electric charge of the struck quark. 



At high momentum transfers, the meson form factor 
can be calculated within the leading order PQCD by 
means of a homogeneous Bethe-Salpeter equation for the 
meson wavefunction. Taking the perturbative kernel of 
the Bethe-Salpeter equation as a part of hard scattering 
amplitude Th, one can get the meson electromagnetic 
form factor given by^ 



F; 



M W ) — 



(fikcPl 

167r3 

.00, 



vI/0?*(y,U)r^(a:,y,qi,ki,li) 



xn^(x,kj_) 

Tc a nV: i±)Th(P{x, k_L), 



(12) 



where Th contains all two-particle irreducible amplitudes 
for J* + qq ^ qq from the iteration of the LFQM wave- 
function with the Bethe-Salpeter kernel. In the 2nd line 
of Eq. H12|l . we combined the spin-orbit wave function 
into the original Th to form a new Th, i.e. 



r - -7? rp{\+X=0) 
^H — K-O-I^H 



rrj „(A + A = ±1) 

'<-±i-Ih ' 



(13) 



where 



TZo = n*^{y,li_)ni;i{x,k±_)+ni^.{y,lj_)ni]ix,kj_) 
_ g [yim2 + y2mi][xim2 + X2mi] 

C'OX Coy 

7^±l = 7^^^(2;,U)7^TT(a;,ki)^-7^J^(2/,U)7^u(x,kJ.) 



Cqx Coy 



(14) 



with Cqx = C and Cqj, = C'{x <-^- y, k ^^ 1^) [see 
below Eq. (O for C]. The hard scattering amplitudes 



T 



(A-|-A=0) 



H 



and T 



(A+A=±l) . 



H 



in Eq. H13|l represent the con- 
tributions from the ordinary-helicity and higher-helicity 
components, respectively. 

To lowest order in perturbation theory, the hard 
scattering amplitude Tnix^y, k±,l±) is calculated from 
the time-ordered one-gluon-exchange diagrams shown in 
Fig-D The internal momenta for (-I-, _L)-components are 
given by 

fci = {xiP+,k^), k2^{x2P+,~k^) 

h = {yiP^,yi<i± + 1_l), h = {y2P^,y2<i± - U)- 

(15) 

In each diagram in Fig. ^ the instantaneous diagrams 
for the intermediate constituents are included using the 
technique shown in Ref. [iJI- In the LF gauge A~^ = 0, 
the gluon propagator is given by 



-9t^y 



k+ 

Kg 



(16) 



^ We should note that the corresponding measure {d? kd? I /IQtt'^] 
in Eq. C3 has to be replaced by [dxd^kx/ievr^lldyd^lx/ieTr^] 
for the BHL-type wave wave function |21| . 




b, Ai D 



D3 A2 D, 



D, A3 D, 






0^2 



Dp 



FIG. 1: Leading order light-front time-ordered diagrams for 
the meson form form factor. 



where 77+ = 0, 77 =1, and ff± = 0. 

Hard scattering amplitudes for the hehcity A + A(=0 
or ±1) components for the diagrams Ai{i = 1,2,3) are 
given by 



T 



(A+A) _ d{y2 - X2) NX 



(A+A) 



Na, 



At 



y2-X2 



DID2 

(A+A) 



T 



(A+A) ^ 6{x2 - y2) N'a > + Na, 
^^ X2 - 2/2 DjDi 

j.(A+A) ^ 0{X2 - V2) A^i'+'' + Na, 
^^ X2 - 2/2 D^D^ 

where the energy denominators are given by 



(17) 



Dx 


Xl 


X2 


D2 


M^ 1 q^ (2/iq±+l-L)^ + m? 
^ 2/1 
(y2q++k+-l+)2 

2/2 -X2 


X2 


D3 


= Di, 




D4 


-M^ + qi (kx+qj^ + mf 

Xl 






(2/2q+ + k+ - l+)2 (j/2q+ - 


l+)2+m2 




2^2 -2/2 


272 


D5 


^^2 k^^ + "^f (2/2q+ + k+ - 
2'i a;2 - 2/2 
(2/2q+ - 1+)^ + ml 


-u)^ 



2/2 
i?6 = i?4. (18) 

The common Na in Eq. (|17|l is obtained from the Feyn- 



man gauge (g'"^) part and given by 



N 



(0) 



2^2 2/1 2/2 qi 



N 



(±1) 



_ -a::ia;2li + 2/i2/2kj 
2;ia;22/i2/2 

-|-X2(xi2/i + 2:22/2)1+ • q+ + 2x22/i2/2k+ • q+ 

-l-(xi2/i + X22/2)k+ • 1+ + xiyiml + X2y2'm\ 

+mim2{yi - Xi){y2 - X2) , 

8[k^ • 1+ + X2I+ • q+ + X2y2ml + xiyiml] 



2^12:22/12/2 



(19) 



where the last mass term 77117712(2/1 — Xi){y2 — X2) in 
N)^ ~ comes from the helicity flip contribution. In 
Eq. lfT7|) . NAi {i = 1,2,3) are obtained from the LF gauge 
parts proportional to 1/fc^ and given by 



Na, 



Na, 



Nas 



2/2 -2:2 



2^2 -2/2 



-{D2 + Di), 
-{D2 + Di), 



X2 - 2/2 



-{D2-Di). 



(20) 



Hard scattering amplitudes for the helicity A 4- A(=0 or 
±1) components for the diagrams Bi{i = 1, 2, 3) are given 

by 



T, 



T, 



(A+A) 
Bi 

(A+A) 
B2 



Ti'+'^(2:^2/,k+^-l+), 



= T 



(A+A) 



T, 



(A+A) _ ^(A+A) 



{x ^ y,k+ ^ -1+), 
\x ^ y,k± ^ -1+). 



(21) 



If one includes the higher twist effects such as intrinsic 
transverse momenta and the quark masses, the LF gauge 
part proportional to l/fc+ leads to a singularity although 
the Feynman gauge part g^^, gives the regular amplitude. 
This is due to the gauge-invariant structure of the ampli- 
tudes. The covariant derivative D^ — d^ + igA^ makes 
both the intrinsic transverse momenta, k+ and 1+, and 
the transverse gauge degree of freedom 5A+ be of the 
same order, indicating the need of the higher Fock state 
contributions to ensure the gauge invariance. However, 
we can show that the sum of six diagrams for the LF 
gauge part(l/fc+ terms) vanishes in the limit that the 
LF energy differences A^, and Aj, go to zero, where A^, 
and Ay are given by 



M^ 



Nr 



ki 



ki 



Xl 



X2 



^ = Af 2 - Ml 



li 



M' 



2/1 



2/2 



Ml. (22) 



In the Appendix, we briefly summarize the proof. 

In this work, we calculate the higher twist effects 
in the limit of A^, = Aj, = to avoid the involve- 
ment of the higher Fock state contributions. Our limit 



Aj; = Ay = 0(but ^/(k^) = f3 y^ 0) may be considered 
as a zeroth order approximation in the expansion of a 
scattering amplitude. That is, the scattering amphtude 
Th may be expanded in terms of LF energy difference 
A as Th = [Tff](o) + A[rff](i) + A^[ThY^^ + • • • , where 
[Th-](°) corresponds to the amphtude in the zeroth order 
of A. This approximation should be distinguished from 
the zero-binding(or peaking) approximation that corre- 
sponds to M = mi + 7712 and k^ — (3 — 0. The point of 



this distinction is to note that [Tf/]*^°) includes the bind- 
ing energy effect(i.e. k_L,lj_ ^ 0) that was neglected in 
the peaking approximation. 

In zeroth order of A^; and Aj^, the net contribution 
from the LF gauge part(l/fc+ terms) vanishes[see Ap- 
pendix] and we only need to compute the Feynman 
gauge(g^'') part, i.e. Na and Nb, for the PQCD analysis 
of meson form factor. The contribution of the Feynman 
gauge to the diagrams Ai is given by 



J 



T 



(A+A) 



1 



Ejn(A+A) ^ j^TJX+X) f S{y2 - X2) 

^^^ ^^ ^ \ y2~X2 D1D2 



N 



(A+A) 



0{X2 - V2) 
X2 -2/2 

A, 



1 



A 



1 



(y2 - X2)DiD2 [ D12 D2 



xy \-i 



A 



(1 



D12' 



(23) 



r 



where D12 = Di — D2 and Axy ~ Ax + Aj,. In the zeroth 
order of A^: and Aj,(i.e. A^. = Aj, = 0), Eq. H23|) reduces 
to 



The hard scattering amplitude for each helicity is sum- 
marized as follows: 



T 



(A+A) 



N 



(A-t-A) 



(2/2 - X2)DiD2 

Similarly for the diagrams Bi, we obtain 



(24) 



T, 



(A+A) 



= T 



(A+A) 



(a; ^ y, k+ ^ -Ij 



(25) 



T. 



(0) 



H 



N 



(0) 



(2/2 - X2)DiD2 



{x <-> 2/, kj 



8 a;2 y 12/2 (2:2 qi + 2k+ • q+) + X2{xiyi + X22;2)l+ • q_L + Xiyim\ 



X2V2m\ + mim2{yi 



xi){y2 -X2) 



2/1 

+ {x 






N 



2/,kj 

(±) 



(a;2qi + 2k+ • c{A_)[x2yl{x2(^^ 



2k+ • q+) - 2x^2/21+ • q+ + (2/2 - X2Ym^] 



+ {x '^ y, k+ 



It 



(2/2 - X2)DiD2 

8 .T2I+ • q+ + X2y2ml H- xiyiml + mim2{yi - Xi){y2 



X2) 



2/1 (a;2qi + 2k+ • q+)[a;22/|(a;2qi + 2k+ • q+) - 2xly2l± ■ q+ + (2/2 - X2)^ml] 



+ (a;^2/,ki^-l+),(26) 



r 



where we neglect the terms such as k^/q^ 



ll/qi, 



and ators and the mass terms (2/2 — X2) 



in the denom- 



can easily find that our result Tjj in the leading twist 
limit reproduces the usual leading twist PQCD result, 
(0) _ iR/^^ „, r^2^ 



]i.± ■ l±/q5^ both in the energy denominators and the nu- inators in Eq. (|26|l give negligible contributions. One 
merators due to the fact that k^^ ^ q^ and 1^^ ^ q'^ 
in large momentum transfer region where PQCD is ap- 
plicable [23 • In the hard scattering amplitudes given 
by Ea. l|2(i|l . the time-ordered 9 function disappears via 
9{x — y) + 0{y ~ x) — 1 and there is no singularity in 
timelike region. We also note that the helicity flip con- 
tributions, i.e. mim2{yi — Xi){y2 — X2), in the numer- 



i.e 



r^' = 16/ix2y2Q'). 



IV. TIMELIKE FORM FACTOR OF A HEAVY 
PSEUSCALAR MESON 

We now consider the PQCD analysis of the timehke 
form factor for the process of e+e" annihilations into 
two pesudoscalar mesons. The hard contribution to the 
timelike form factor for the electron-positron annihila- 
tions into two pseudoscalar mesons, i.e. e+e~ — > MM, 
is obtained as 

Fniq^) = ei/(g^,TOi,m2) -f 62/(5^,7712, mi), (27) 

with the amplitude I{q^,mi,m2) given by 



powers of (kj_ • q_L) and (1_l • q±) in the numerator since 
TZo[TZ±] includes even[odd] powers of k_L and 1_l. As a 
result, we can get Th in Eq. 1)28(1 as a function of even 
powers of k^ and Ij^ for both ordinary and higher helic- 
ity contributions. We then analytically continue to the 
timelike region by changing q^ to iqj^(or q^ -^ —q^ in 
this case) in the form factor. 

In terms of FM{q^) given by Eq. (|27|l . the cross section 
of the pseudoscalar meson pair(MM) production in the 
unpolarized e"'"e~ annihilations is given by 



da 
dh 



(e+e" -^ MM) 



3^ 
327r 



2m2 



<T^+e-_+ - sin- 0\FM{q )| 



/(9^ 



mi, 777,2 j 



nCpas 



dxdyd'^k±d'^l± 



167r3 

x0(y, l±)XH0(a;,k_L), 




(30) 



where /? 



^1 - 4M2/g2 and (t^+. 



^n-i-fj.- 



where as is the QCD running coupling constant and 
Cp{— 4/3) is the color factor. We note that all the in- 
variant masses in TZq and TZ±i in the spin-orbit wave 
function are replaced by the physical meson mass to be 
self-consistent with the result for the hard scattering am- 
plitude in zeroth order of the LF energy differences A^, 



^«V(3i?Lam) with Sbeam = vV/2. 

In the peaking approximation, the transverse momenta 
of the quark and antiquark are neglected and the longitu- 
dinal momentum fractions are given by xi = yi = mi/M 
and X2 = y2 — m^jM with M = mi + 7712- In this ap- 
proximation, the higher helicity contribution to the hard 
scattering amplitude also vanishes and the ordinary he- 
licity contribution to the hard scattering amplitude given 
by Eq. ((26|l can be rewritten as 



and Ay, i.e 



^x ) ^y 



0. Again, our analysis should 
be clearly distinguished from the peaking approximation 
(zero-binding limit or /3 ^ 0) which leads to the DA 
(J)m{x) defined by Eq. lO as the (5- type function with 
Xi = rrii/M. In our zeroth order approximation of LF 
energy differences, we consider /3 7^ (i.e. the effect 
of the transverse size of the meson). In the next sec- 
tion (Sec.V), we take the /3 values determined from the 
analysis of the mass spectroscopy using our LFQM with 
the variational principle for the QCD-motivated effective 

>.^ , 1^ ciiiui n.^ - ij_ terms compared 



i^H J peaking 



-xlviy^q^ + xiyiml + X2y2m\ 



xlvivl 

+ {x ^ y) 
32M7 / M 
g4 \^TO2 



1- 



q^ 2m2 



4M2 mi 



(31) 



Hamiltonian 



alpr 



where 7 = 777 1 7772 /Af. 

Therefore, the peaking approximation of the timelike 
form factor of a heavy pseudoscalar meson is given by 



Since we neglect the k^ , 1^ and k 



to large q^ as we stated below in Eq. ^^, we have only [/Af]pcaking(9^) oc ei / dxdyS{xi - mi/M)[Tjj ]pcak 



k_L(l_L) ■ q± and qj^ in both numerator and denominator 
terms in Eq. (|26|l . Thus, for convenience in our numer- 
ical calculation, we change the denominator in Eq. I|26|) 
to include only even powers of k^ and 1_l. Then the 
numerators of the hard scattering amplitudes T^ and 
T^ ' in Eq. ((^ include both even and odd powers of k_L 
and lj_ via the terms (kj_ • qj.)™ and (1_l • qx)™, where 
777 =integer(non- negative). After the change, the generic 
form of the hard scattering amplitude may be given by 



xS{y^-m,/M) + (1^2) 



oc ^iei 



+ 62 



M 

mi 



Af ^ 


4 
\ 


m2y/ 


^ 4 






1- 


y ) 





2m2 



4A/2 m,i 



q^ 2mi 



4M2 



m2 



(32) 



Th^ 



iV((k^ • q^)", (U • q ^)'") 
i?(qi",(kx-q±)2«,(lx 



This reproduces the result obtained by Brodsky and 
Ji 11]. The form factor zero in this approximation occurs 
at 



\-2n\ 



(29) 



where 77 and ?7i are (non- negative) integers. In Eq. ((29|l . 
we show only essential terms, kj^ • qj^ and lj_ ■ q_L in 
the numerator to explain how to obtain the nonvanishing 
contributions from the ordinary and the higher helicity 
contributions. That is, as one can see from Eqs. H13|) 

and IHI), by combining T^°^[T^='='] in Eq. ^ with the 
spin-orbit wave function 7^o[''^±] in Eq. (fT^ to get Th 
in Eq. USJ or Eq. (gHl, t''h\t''h'^] should have even[odd] 



mi 
2m 2 



ml 
2mY 



4Af2 



1+ ^ 



(33) 



Even though 7711 > 7772, the 62 contribution is not negligi- 
ble for the heavy-heavy pseudoscalar meson system such 
as Be- The reason for this is because the timelike form 
factor of heavy pseudoscalar meson encounters a zero and 
62 contribution in the region near the form factor zero 



has non-negligible effect. However, for the heavy-light 
quark system such as B and D mesons, the light quark 
contribution(i.e. 62) can be safely neglected. 



V. NUMERICAL RESULTS 

In our numerical calculations, we use the model pa- 
rameters (rnQQ, Pqq) obtained from the meson spec- 
troscopy with the variational principle in our LFQM [ig | 
for the linear confining potential. Our model parame- 
ters are summarized in Table I. As mentioned earlier, we 
should note that the r.m.s value of the transverse mo- 
mentum in our LFQM is equal to the gaussian (3 value, 

i-e- y^(ki)QQ=/3QQ. 

The shape of the quark DA which depends on the /? 
value is important to the calculation of the cross section 
for the heavy meson pair production in e+e~ annihila- 
tions. We show in Fig. |5] the normalized quark DA of 
B and Ds[D] mesons with different values of (3. For 
the quark DA of i3-meson, we compare our LFQM re- 
sult (solid line) with the small f3 value result close to the 
peaking approximation, e.g. /3 = 0.1 GeV(dashed line). 
As one can see, our LFQM result for the quark DA of B 
meson shows sizable deviations from the peaking approx- 
imation. For the quark DA of D^ [solid line] and £> [dotted 
line] mesons, the peak for (J)d{x) is located to the right of 
that for (j)r)^{x). This indicates that the c-quark carries 
more longitudinal momentum fraction in D than in Ds 
as one may expect. We also show the ratio [dash-dotted 
line] of (l>D^{x) and (J)d{x) for the sake of comparison. 

In Fig. O we present also the normalized DA of vari- 
ous ^5*0 quarkonium((7g) states obtained from our LFQM 
parameters given in Table I. To explain the discrep- 
ancy between the NRQCD prediction [3] and the Belle 
measurement B for (T(e+e~ -^ J/il^ + Vc), Bondar and 
Chernyak(BC) |Tj] reasoned that the discrepancy may 
be due to the extreme (5-function-like charmonium DA 
adopted from NRQCD and claimed that they can fit the 
Belle data by choosing a rather broad DA for the char- 
monium state. Interestingly, our LFQM prediction for 
(pec shown in Fig. O looks quite similar to BC's result 
in a sense that the DAs for heavy quarkonium states 
differ from the (5-function-like DA. In our model calcula- 
tion, the DA gets narrower as /3 gets smaller. Also the 
timelike form factor Fmiq^) with small (3 value decreases 
faster than that with large /3 value. Since the cross sec- 
tion (Te+e-^jvfM i^ proportioual to |Fm(9^)P, the cross 
section with small f3 is small compare to that with large 
/3. Thus, the cross section for e+e^ -^ MM can be in 
principle enhanced by broadening the quark DA. 

In Fig. ^ we compare the results of (t)cc{x) in more 
detail. The solid and dashed(dash-dotted) lines repre- 
sent our LFQM result and BC's result Q with w^ = 
0.15(u^ = 0.3), respectively, where the parameter v\l^ 
represents the characteristic quark velocity in the bound 
state. Although there is a similarity in the quark DA be- 
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FIG. 2: Normalized quark distribution amplitudes of _B[top] 
and Ds [D\ [bottom] mesons for different values of the gaussian 
parameter (5. 



tween ours and BC's results(in particular, their u^ = 0.15 
result), there is a rather substantial difference near the 
end-point region between ours and BC's results. Since 
the PQCD hard scattering amplitude is typically very 
sensitive to the end-point values of DA, it may not be so 
difficult to imagine that BC's prediction of double charm 
production cross sections would have made a difference 



TABLE I; The constituent quark masses mg[GeV] and the gaussian paramters l3{- 
obtained from the variational principle. q=u and d. 



^/{kjy)[GeV] for the linear potential 
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FIG. 3: Normalized quark distribution amplitudes of qq with 
our LFQM parameters f3qq. 



FIG. 4: Normalized quark distribution amplitude of cc(solid 
line) compare to those(dashed, dash-dotted lines) obtained 
from Bondar and Chernyak [l^. 



depending on what w^ value they have used. To fit the 
Belle measurement yf for a{e~^e~ —> J/tp+ilc)! they used 
w^ = 0.3 rather than v'^ = 0.15. However, as we stressed 
earlier, the way that BC[lJI handled the transverse mo- 
mentum effect is different from ours since they integrated 
out the transverse momentum seperately in the wave 
function part and in the hard scattering part while we 
didn't factorize the hard and soft parts but integrated out 
the transverse momentum for the whole amplitude. We 
think that a consistent analysis with /3 ^ (or non-zero 
binding energy) should follow our non-factorized formu- 
lation (see e.g. Ea. lf^ ). 

In Fig. [3 we show the cross sections for the exclu- 
sive DfDj and D^D~ pair productions in e+e~ anni- 
hilations. The solid line represents the results including 
both ordinary and higher helicity contributions while the 
dotted line corresponds to the result of the ordinary he- 
licity contribution only. The dashed line for (T(e^e~ — > 
D+D") represents the lower limit(i.e. (3 = 0.16 GeV) 
for the form factor zero to occur above -y/s = 10.6 
GeV. The small black circle for (T(e+e^ -^ D^D^) 
represents the upper limit obtained from Belle |2l, i.e. 
o-Exp.(e+e- -^ D+D-) < 0.04. The higher helicity 
contribution to the cross section, i.e. the difference 
between solid and dotted line, is more pronounced in 



e+e~ -^ D^D^ than in e+e^ -^ DfD^ especially 
near the turnover point (or the form factor zero point). 
In general, the higher helicity contribution to the me- 
son form factor increases(decreases) as quark mass de- 
creases (increases). For instance, while the higher helicity 
contribution to the hard scattering amplitude is negligi- 
ble in the B meson form factor, it is not negligible in the 
pion form factor. However, the most significant in our 
analysis is the transverse momentum effect which delays 
the turnover point (see Eq. H33|l for the peaking approx- 
imation). For instance, the turnover for Ds[D] meson 
occurs near 6. 3 [11. 3] by going beyond the peaking ap- 
proximation while the corresponding turnover point is 
near ^jM ~ 2.8 [4.03] for the peaking approximation. 

Numerically, we obtain the cross sections for D^D^ 
and D~^ D~ pair productions at ^/s = 10.6 GeV with our 
LFQM paramters(i.e. /3 = 0.5016 GeV for Ds and 0.4679 
GeV for D) as follows 



a(e^e 



DTDs 
D+D- 



{S.Otti) X 10-4[pb], 
(0.02 ± 0.01) [pb]. 



(34) 



for the strong coupling constant a, 
lar ttc values were used in Refs. |2 



0.2 ±0.05. Simi- 
12 J . Our result for 



U) 


K' 


, , , 1 , , , , 1 , , 


= 




: \ 




_ *LPOM(''.M;P=0-5016.Gey 
(oral nary +nig her hehcities) 


: 


10-' 






fLFOM(''.V,P=0-5016GeV 
(ordinary nehcity only) 


= 




l'2 , r 


.6GeV 






10'^ 


E" A 




— 


















































a 


\ 




" 


D 10-' 


— --A 




^ 




= 




: 


*v 


- 


^ 


- 




" 


\\ ^_;,_,,^__^ 


" 


10"* 


r 




1 


10-' 


E- 




-s 


.^-6 


1 




^ 



13=0.4679 GeV 
|3 = 0.16GeV 
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FIG. 5: Cross sections for DfD^ [top] and D^D [bottom] 
pair productions in e^e~ annihilations with Os — 0.2. 



a{e~^e~' — > D^D^) is consistent with the recent experi- 
mental data from Belle, (TExp.(e^e^ -^ D^D~) < 0.04. 
Since a{e'^e~ -^ D^D~) gets larger as /3 grows, the up- 
per bound of (TExp. from Belle provides a constraint on 
the maximum /3 value modulo the dependence on as- 
If the cross section for the D meson pair produc- 



tion satisfies a{e^e 



D+D-) < 0.04 and its slope 



with respect to the momentum transfer is negative, i.e. 



{da/dy/s) < at \/s — 10.6 GeV, then we could also set 
the lower bound for /3 value as /3 > 0.16. The shape of 
the D meson quark DA corresponding to /3 = 0.16 GeV 
is shown by dashed line in Fig.|B| If /? < 0.16 GeV, the D 
meson quark DA approaches to the J-type function but 
(dcr/d^/s) > at Y^ — 10.6 GeV due to the form factor 
zero occuring at ^/s < 10.6 GeV. Due to the occurence of 
form factor zero for the heavy pseudoscalar meson pair 
production |j, IllL l25l |. more experimental data around 
^/s = 10.6 GeV are necessary to check the slope of the 
cross section. More data around ^/s = 10.6 GeV would 
further constrain the shape of D meson quark DA. 

How about B mesons? We should point out that the 
PQCD result of the cross section for B^B^ pair produc- 
tion may not be trustworthy because y/s = 10.6 GeV is 
too close to the threshold energy oi B'^B~ pair produc- 
tion. As expected, the gluon momentum transfer from 
the heavy quark to the light quark in B'^ B~ pair produc- 
tion at ^/s — 10.6 GeV turns out to be only around a few 
hundred MeV close to the scale of Aqcd ■ On the other 
hand, the gluon momentum transfer in DT^, £),~, pair pro- 
duction at y^ = 10.6 GeV is much larger than the scale 
of Aqcd ■ By going beyond the peaking approximation, 
the average gluon momentum transfer gets even larger 
due to the transverse momentum effect. This may jus- 
tify our PQCD analysis for _Dr^,_Dr, pair production at 

V^ = 10.6 GeV. 

Although the absolute value of the cross section for 
i3[i3s]-nieson may not be reliable near ^/s = 10.6 GeV, 
it seems interesting to discuss the behavior of the ratios 
of cross sections such as a{e~^e~ -^ B^ B'^) / a{e^ e~ — > 
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FIG. 7: Our predictions [closed symbols] on the cross section 
ratios for various heavy pseudoscalar meson(i?, Bs, -D,and Ds) 
pair productions compared to the peaking approximation re- 
sults[open symbols] near ^/s — 10.6 GeV. 



B+B^) and a{e+e- 



B,"B,")/a(e+e- 



B+B- 



In Fig. 13 we show our predictions [closed symbols] on 
the cross section ratios for various heavy pseudoscalar 
nieson(i?, i3s, Djand Ds) pair productions in e+e~ an- 
nihilations near -^s = 10.6 GeV, i.e.^/s ranging from 
10 to 12 GeV and compare our results with those [open 
symbols] obtained from the peaking approximation [ill ]. 
The open and closed diamond symbols for cr(e^e~ — > 
B^B^)/a{e^e~ -^ B^B^) are on top of each other and 
their values are almost equal to 1. In fact, the cross sec- 
tion ratios involving light quarks u and d such as (-D*^, D) 
and (-B*^, B) cases are close to 1. This is due to the negli- 
gible contribution from the diagrams where the photon is 
attached to the light quarks. As the replacement of light 
quarks by the strange quark makes those diagrams non- 
negligible, the cross section ratios for the cases of {Ds^ D) 
and (Bs,B) deviate from 1 appreciably. However, most 
significant is again the transverse momentum effect which 
is pronounced in the case of D[Ds] meson pair produc- 
tions compared to B[Bs] meson pair productions. In par- 
ticular, the deviation between the open and closed sym- 
bols for the case {Dg , D) is quite dramatic compared to 
the case of {B^.B). 



found to be relat ed to the transverse momentum via 
Pqq — J (k\)QQ. This relation naturally explains the 

zero-binding energy limit for the zero transverse momen- 
tum, i.e. (Mq) — {mi + 7712)^ and Xi — rrii/M for (3 — Q. 

However, the heavy quark DA is sensitive to the value 
of /3 and indeed substantially broad and quite different 
from the (5-type DA according to our LFQM based on 
the variational principle for the QCD-motivated Hamil- 
tonian |16l \u\ . If the quark DA is not an exact 5 func- 
tion, i.e. k_L in the soft bound state LP wave function 
can play a significant role, the factorization theorem is 
no longer applicable. To go beyond the peaking approx- 
imation, the invariant amplitude should be expressed in 
terms of the LF wave function ^(x^, k^,;) rather than the 
quark DA. 

In going beyond the peaking approximation, we 
stressed a consistency by keeping the transverse momen- 
tum kj^ both in the wavefunction part and the hard scat- 
tering part together before doing any integration in the 
amplitude. Such non-factorized analysis should be distin- 
guished from the factorized analysis where the transverse 
momenta are seperately integrated out in the wavefunc- 
tion part and in the hard scattering part. Even if the used 
LF wavefunctions lead to the similar shapes of DAs, pre- 
dictions for the cross sections of heavy meson productions 
would apparently be different between the factorized and 
non-factorized analyses. 

In this work, we compared our non-factorized analy- 
sis with the usual factorized analysis based on the peak- 
ing approximation and found a substantial difference be- 
tween the two in the calculation of the cross section for 
the heavy meson pair production. We also discussed the 
higher helicity contribution to the cross section. Our 
analysis provided a constraint on the size of quark trans- 
verse momentum inside the D meson from the recent 
Belle data, o-Exp.(e+e" -> D+D') < 0.04. More ex- 
perimental data around y/s — 10.6 GeV would further 
constrain the shape of D meson quark DA and test 
our LFQM prediction. Application of our non-factorized 
PQCD analysis to the alphas higher order corrections, 
e.g. in double charm production, would deserve further 
investigation. 
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APPENDIX A: PROOF OF VANISHING 
LIGHT-FRONT GAUGE PART IN M = Mo LIMIT 



obtained from Eqs. H17|) and (|20|l is given by 



The contribution of the LF gauge part(l//Cjj' term) to 
the hard scattering anipHtude for diagrams A ~ X]i=i -^i 



T 



(i/fe+) S{V2 - X2) Na^ , 0{X2 - V2) ( Na 



V2 --X2 



D1D2 
D2 + 



X2 -V2 



D, 



(2/2-2:2)2 D1D2 



Na, 



D^D. D.D, 



d{V2 - X2) +9{X2 - 2/2) 



r 



D2-DiDiD2 + D^{D2 + D^) 



D4 



{D2 + Di)D5 



(Al) 



In terms of the LF energy differences, the relevant energy 
denominators can be rewritten as 



Di = 



D2 + Di 
D2 + D5 



1 / 2 

Xl 

A,, 



2ki • q±) 






A„ 



(A2) 



Eq. IjAip and the corresponding LF gauge part to the 
diagrams Bi lead to singularties. Fortunately, however, 
in the zeroth order of A^, and Aj^, i.e. A^;, Aj, — > limit, 
one can see that the energy denominator term in 9{x2 — 
2/2) in Eq. (|A1|) vanishes, which leads to 9{y2 ~ X2) + 
0(x2 — 2/2) — 1- Thus, the LF gauge part contribution to 
the diagrams A — X]i=i ^i becomes 



T 



(i/fc^) 



(2/2 - X2Y' D2 



(A3) 



and similarly we obtain 

^(i/fc+) 



(2/2 
B, 



X2Y Ds 



(A4) 



for the diagrams B = X)i=i ^i ■ Finally, from the rela- 
tion D2 + Dg ~ Aa; + Aj^, one can see that the net contri- 



bution from the LF gauge parts, i.e., T 
vanishes exactly in the limit of Aj, 



0. 
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provides a constraint on the size of quark transverse momentum inside the D meson from the recent 
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Belle data, (TExp.(e/-TS_.^,i?^J2. 
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[= peaking approx. by Brodsky and Ji 



I. INTRODUCTION 



1 




the rtheoretical predictior»i-|6E)b(TEti ^</^5f^ T|^8]\ and 
)cesses although higher orocr correcx: 



incl isive [9] processes alttioilgti tiigtiet order corrections 

-reduce the differences [10]. A couple of years ago. 

Belle Collabotation [3] reported the first measure- 

t of the e+e" -> D*+D*~ and e+e~ -^ D+D*~ 

sections and polarizations ^t ^^ — 10-6 GeV. ^ 

Thc|y also set an upper lj>eytT>n fee civn^^(BtidE> f^r 

D^D^ . Interestingly, wkile tfee T;neorcEicaTprc- ^ 
ions based on the heavy quark effective theory [4] and 
onstituent quark model [7] for e+e" — > D*'^D*~ and 
-^ D^D*^ cross sections are similar to the mea- 



l(jfeurc3~data [3], the predictions for e"* 



cross 



sect ton are either, 
than the dat 



ither nuitte smaller [4j or sfemfiwl 



D+D- 

hat larger [7] 



Tfte aWve e^lus^e/inclusive meson pair productions 
pro^ade a unique O]^orti0ity ^ investigate asymptotic 
behiviors of various meson form fJfcorswn the frame- 
wor ~of perturbative quantum chromodynamic^fPQflp). 
The heavy meson pair production is of special interest ' 
sinc3"gluons carrying large momentum transfers can be 
rati er easily accessible in the kinematic region above the 
threshold. Also, the wavefunctions of heavy systems may 
b2 ^1^1 constrained due to the heaviness of constituents. 
IL/This-, it has been pointed out that exclusive pair pro- 
duc Ton of heavy mesons can be reliably predicted within 
PQCD [11]. 

Iflhc factorization theorem in PQCD is applicable to 
exclttsive processes, then the invariant amplitude for ex- 
clusive processes factorizes into the convolution of the 
valence quark distribution amplitude(DA) (l){x,q^) with 
the liard scattering amplitude Th [12]. To implement the 
fact Drization theorem at high momentum transfer, the 
had'onic wave function plays an important role linking 
bet^reen long distance nonperturbative QCD and short 
d3;tknco PQCTD. A f)a,rticliila,rl4 con-^^enieilt andl intniltive 



10 



framework in applying PQCD to exclusive proces'Sjs is 
based upon the light-front (LF) Fock-state deconi])osi- 
tion of hadronic state. In the LF framework, the va- 
lence quark DA is computed from the valence LFji^ave 
function \['„(a;i, k^j) of the hadron at equal LF :ime 
t + z/c which is the probabilibity amplitude to find 
t)nstJBient^quaR;s,aiHiqu£HiS, Hid ^lon^ wiiS LF 
_ snt^j =^a;,^,)^a^dr^ Ej^, .^aH^^_Lj 
arc the LF momentum fraction and the transverse mo- 
lenta of the ith constituent in the n-particle Fock--s ;ate, 
^spectively. a A A A / 

?o lowest o^dcr Ak powurbarion ilieory of the m 3son 
fo^ faCTor calculation at large momentum transfei^ the 
hard scattering amnlitud^ ThJ§ doj^inalAl b5wne-4 luon 
exchan^dia^ams. For the factorization theorem Id be 
apP.icable in the heavy meson pair production analj^sis, 
the only consistent form of the quark DA would be t he 5 
function, i.e. (j){x, q^) ~ S{x — uiq/M) where uiq an 1 M 
are the heavy quark mass and the meson mass, ret pec- 
tively [13]. In this so called "peaking approximation", 
the momentum fraction carried out by ith constjl lent 
is equal to the ratio of the constituent mass to mjson 
mass, Xi = TTii/M. This relation implies M = mi T m2, 
i.e. zero-binding energy limit. ~ 

However, as pointed out in Ref. [13], if the quark DA 
is not an exact S function, i.e. kj^ in the soft be und 
state LF wave function can play a significant role, the 
faa4prization theorem is no longer applicable. To ge be- 
yond thl^)eaking approximation, the invariant ampli ;ude 
should be ex™ess^ in terms of the LF wave fune tion 
^(a;i,k^i) rather than ^e quark DA. In Ref. [13], the 
validity issue of peaking apprc»jmation for the heai^y me- 
son pair production processes was discussed using tE ; LF 
model wave function 'i'{xi,'k±i) oc *p(-^J^o//3^), where 
Mo is the invariant mass of the constituent quark and 
antiquark defined by Mq = X^iC^l + rnf)/3^an(S' /3 is 
the gaussian parameter. The limit /3 — > corresponds 
to the peaking approximation(i.e. zero-binding eyrgy 
limit M = mi + 7712). In the analysis of the heavy-h javy 
system like Bc{bc) meson, it was found that the effect of 
going beyond the peaking approximation (/3 up to 100 
MeV) was not important compared to the peaking ap- 
proximation limit(i.e. /3 r-> 0) [13]. However, it is no ; yet 
clehT if the s^me chnchiMnn -iKfouldl applt/ to the hiUa.vy- 
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10.5 



11 

1/2 



11.5 



12 



light system such as D and B mesons. Moreover, the 
initial analysis limited only up to /3 < 100 MeV may not 
be sufficient to draw a definite conclusion on the validity 
of the peaking approximation. 

The main purpose of this work is to extend the previous 
analysis [13] and point out that the recent Belle data [3] 
can provide a rather stringent constraint on how broad 
or narrow the D^ meson quark DA is. Clarifying the 
relation between the /3 value and transverse momentum, 
kj^ , is a particularly important issue since the quark DA 
is very sensitive to the (3 value and the different shape of 
the quark DA could enhance or reduce the cross section 
for the exlcusive meson pair production in e^e~ annihila- 
tions. Incidentally, Bondar and Chernyak [14] considered 
a rather broad quark DA(i.e. rather significant binding 
energy effect) instead of (5-type quark DA to explain the 
data for the exclusive e+e" -^ J /ip + rjc process. Ma and 
Si [15] also previously discussed the variation of DA to 
explain the data for the same process. 

In this work, we stress a consistency of our analysis 
in going beyond the peaking approximation. In particu- 
lar, we confirm that the (3 value in our model LF wave 
function is related with the transverse momentum via 
(^QQ — Vv^^jTqq- ^^ expected, the non-zero (i value cor- 
responds to the transverse size of the meson and /3 — > 
limit corresponds to the peaking approximation(i.e. zero 
binding energy limit) as discussed in [13]. This implies 
that it may be significant to keep the transverse momen- 
tum k_L both in the wavefunction part and the hard scat- 
tering part together before doing any integration in the 
amplitude if (3 is not so close to zero or the binding en- 
ergy effect is not negligible. Thus, we think that the 
factorization of amplitude by integrating out the trans- 
verse momentum seperately in the wavefunction part and 
in the hard scattering part cannot provide a consistent 
analysis to take into account the binding energy effect. 
This may distinguish our method from Ref.[14] to take 
into account the binding energy effect. 

We also note that our gaussian parameter (3 is not 
chosen arbitrarily but fixed by the variational principle 
for the well-known linear plus Coulomb interaction mo- 
tivated by QCD [16], which in turn uniquely determine 
the shape of the quark DA in our model calculation. This 
implies that the recent data by the Belle collaboration [3] 
provide a useful test on our model calculation. 

The paper is organized as follows. In Sec. II, we 
describe the formulation of our light-front quark model 
(LFQM), which has been quite successful in describing 
the static and non-static properties of the low- lying 
mesons [16, 17]. In Sec. Ill, the transverse momentum 
dependent hard scattering amplitude for the meson is 
given within the LF framework. The contribution to the 
meson form factor from higher helicity components is also 
given in this section. In Sec. IV, the analytic continuation 
from spacelike region to the timelike region is introduced 
to obtain the cross section for the pseudoscalar meson 
pair(MM) production in e+e^ annihilations. As a valid- 
ity check of our model, we also show that our result for 



the meson form factor obtained in Sec. Ill reduces to the 
peaking approximation in the (3 — > 0(i.e. zero-binding) 
limit. In Sec. V, we present the numerical results for 
the e+e^ -^ D'^,D7, cross section and compare with the 
available data. Summary and conclusions follow in Sec. 
VI. In the Appendix, we briefiy summarize our proof of 
vanishing contribution from the light-front gauge part in 
M = Mo limit. 



II. MODEL DESCRIPTION 



In our LFQM, the meson wave function is given by 



*5[/^(x,ki,AA) 



{x,k^)nl'l^ix,k^). (1) 



where 



fa;,ki 



is the radial wave function and 



TZ.T'{x,'k±) is the spin-orbit wave function obtained by 
the interaction-independent Melosh transformation [18] 
from the ordinary equal-time static spin-orbit wave func- 
tion assigned by the quantum numbers J^*^ . The meson 
wave function in Eq. (1) is represented by the Lorentz- 
invariant variables Xi ~ p'l /P^, k^j = p±i — XiP± and 
Xi, where P,Pi and Xi are the meson momentum, the 
momentum and the helicity of the constituent quarks, 
respectively. 

The radial wave function (f>(x,]i.±) of a ground state 
pseudoscalar meson (J-^*^ = ^) is given by 



(x.ki 



1 



7r3/2/33 



1/2 



exp(-fcV2/3'), (2) 



where the gaussian parameter /3 is related with the size 
of the meson. Here, the longitudinal component k^ of the 
three momentum is given by k^ = {xi — ■k)Mo + {ra\ — 
ra\)/2Mi) with the invariant mass 



M,i 



ki 



ki 



ml 



(3) 

Xi X2 

where Xi = x and a;2 = 1 — a;. The covariant form of the 
spin-orbit wave function 7?.°^(x, kj^) for the pseudoscalar 
meson is given by 






u{pi,X)lbv{p2,X) 



V2[Mo2-(mi-m2)2]i/2' 
and its explicit matrix form is given by 



(4) 



.^00 ^ _L 



-k^ 

-Xim2 — X2mi 



XiJn2+X2mi . , , 



where C = ^y2xlX2[MQ - (mi - 1712)^] and k^^'^^ 



ik- 



y. Note that Eaa<>a"a = 
our wave function is given by 



rCry. zh 



1. The normalization of 



^ fd^k\^"j^{x,k^,XX)f 

AA •' 

= J\xjd'k^(^^ycbix,k^r^l, (6) 



where the Jacobian of the variable transformation 
{x, k^} ^ fc = (kj_, fc^) is given by 



dx 



Mq 

4xia:;2 



1- 



(mi - 7712) 



212 



M§ 



(7) 



The effect of the Jacobi factor has been analyzed in 
Ref. [19]. 

With this normalization, the root-mean-square(r.m.s.) 

value of the transverse momentum, i/(k^)gQ is obtained 

via 

(ki)QQ - eg [d'k\klMx,k^)\^ + {Q^Q). (8) 



Numerically, we confirm that w(k^)QQ = f^nQ. The 
numerical values of (3qq are discussed in Sec. V (see 
Table I). We also find the relation {M§) > (mi + 7^2)^ 
for (3 > 0. Our numerical solution indicates that the 
zero-binding energy limit, M = Mq = mi + m2, is ac- 
complished by taking both k^ ^ and /3 ^ limits in 
our LFQM. 

The quark distribution amplitude(DA) of a meson, 
<l>M,x{x,Q), i.e. the probability of finding colHnear quarks 
up to the scale Q in the L^ = O(s-wave) projection of the 
meson wave function [12] is defined by 



't'M,\ix,Q) 



Q 



[d2kJ*(.T,ki,AA), 



(9) 



where [d^k^] ^ (Pkj_^/dkJd^/Vm^ for * = *5?. 



III. HARD SCATTERING AMPLITUDE WITH 
kx-DEPENDENCE 

In this section, we calculate the pseudosclar meson 
electromagnetic form factor in the region where the 
PQCD is applicable. Our calculation is carried out using 



the Drell-Yan-West frame [20] (g"" 



0) with 



q^ — Q^ — —Q^- The momentum assignment in the 
g+ = frame is given by 

P = (P+,^,0J, P' = (P+,^^^-±^,q^) 



where d^k — dxd^k^^Jdk^/dx, k'^ = k^ +X2q± and Cj 
is the electric charge of the struck quark. 

At high momentum transfers, the meson form factor 
can be calculated within the leading order PQCD by 
means of a homogeneous Bethe-Salpeter equation for the 
meson wavefunction. Taking the perturbative kernel of 
the Bethe-Salpeter equation as a part of hard scattering 
amplitude Th, one can get the meson electromagnetic 
form factor given by^ 



T;iHard^,o2N 

-T J, 



M 



'{Q' 



d^kdH 



167r3 

,00/ 



vl/r(y,l±)rH(x,y,q^,ki,li) 



xM/^(a;,kx) 

Qj rZQ, t , Xz-r / / 1 \ 

,^ , (P\y, 1±)Th(P{x, k_L), 
Idtt-^ 



(12) 



where Th contains all two-particle irreducible amplitudes 
for J* + qq -^ qq from the iteration of the LFQM wave- 
function with the Bethe-Salpeter kernel. In the 2nd line 
of Eq. (12), we combined the spin-orbit wave function 
into the original Th to form a new Th, i.e. 



T - -7? rpiX+X=0) 



T-, „(A+A=±1) 
l<-±lJ^H ' 



(13) 



where 



7^o = 7^^x(y,U)7^Ti(x,ki)+7^J1-(y,U)7^xT(a;,ki) 
_ g [Vim2 + V2mi\[xim2 + X2mi] 

7^±l ^ 7^^^(y,U)7^Tt(x,ki)+7^;[^(y,U)7^u(2;,ki) 



Cqx Coy ' 



(14) 



with Cqx = C and Coy = C{x ^^ y,k ^^ Ij^) [see below 

Eq. (5) for C]. The hard scattering amplitudes T^ ~ ' 

and T^ ~ in Eq. (13) represent the contributions 
from the ordinary-helicity and the higher-helicity com- 
ponents, respectively. 

To lowest order in perturbation theory, the hard 
scattering amplitude rfr(a;,y, k^, Ij^) is calculated from 
the time-ordered one-gluon-exchange diagrams shown in 
Fig. 1. The internal momenta for (-f ,_L)-components are 
given by 



Q = (0,pt,q±), 



(10) 



where prime denotes the final state momentum and q = 
P' — P and M is the physical meson mass. 

As a starting point, the electromagnetic form factor of 
a pseudoscalar meson is given by a convolution of initial 
and final meson wavefunctions: 



fcl = {xiP+,k^), k2 = ix2P^,-kj_) 

h = (yi-Pi^, 2/iq± + U), h = {v2Pi,y2C{i_ - U). 

(15) 

In each diagram in Fig. 1, the instantaneous diagrams 
for the intermediate constituents are included using the 



FTm = ^5]e,y'd3fcvI>°o*(x,k'^,AA) 



xM'™J(x,ki,AA), 



(11) 



^ We should note that the corresponding measure [(fikd^l/Wn^] 
in Eq. (12) has to be replaced by [dxd^l(.j_/16n^][dyd^lj_/16n^] 
for the BHL-type wave wave function [21]. 






D '^' D 



FIG. 1: Leading order light-front time-ordered diagrams for 
the meson form form factor. 



The common Na in Eq. (17) is obtained from the Fcyn- 
man gauge((7^'') part and given by 



N 



N 



(0) ^ 22 i2 1 2 

xiX2yiy2 I 

+X2ixiyi + a:2j/2)l± • q± + 2x2yiy2^± ■ q± 

+ {xiyi + X22/2)kj_ • 1± + xiyiml + X22/2™? 

+Tnim2{yi - xi){y2 - X2) , 

(±1) _ 8[k^ ■ U + a:2l± ■ q_L + X2y2ml + xit/ito^] 



a;ia;2yi2/2 



(19) 



where the last mass term mi7Ti2(yi — xi){y2 — X2) in 
Ny^ ~ comes from the helicity flip contribtion. In 
Eq. (17), NAi {i = 1,2,3) are obtained from the LF gauge 
parts proportional to 1/fc^ and given by 



technique shown in Ref. [12]. In the LF gauge A~^ = 0, 
the gluon propagator is given by 



^^i^ y^i^ I 



{kg)^,r|^ + {kg)i,T]^, 



kg 



(16) 



where 77+ = 1, r;^ = 0, and fj±^ ~ 0. 

Hard scattering amplitudes for the helicity A + A(=0 
or ±1) components for the diagrams Ai{i = 1,2,3) are 
given by 



(A+A) 



T 



(A+A) _ 0{y2 - X2) JVr ^ + Na, 



y2 - X2 D1D2 

r(A+A) 



T 



(A+A) _ 0{x2 - y2) A^r ^ + Na, 



X2 - 2/2 D^D4, 



T 



(A+A) ^ 6{x2 - y2) N'r > + Na, 

"^'^ X2 - 2/2 D^^Dq 

where the energy denominators are given by 

^^^^2^q2^_(k++q+)^ + m? ki 



(17) 



D2 = M2 + q2 



Xi X2 

(?yiq++U)2 + m2 k^ + mi 



2/1 



X2 



(2/2q+ + k+ - 1+) 



2/2 - X2 



D3 - Di 

D4 ^ M' + qi- 



2 , „2 (k+-fq+)2 + m? 



a;i 



(2/2q±+k+-l+)2 (y2q+-l+)2-L-2 



7T15 



2:2 -2/2 



2/2 



D, ^ M^ - 



2 k\+ml (y2q+ +k+ -l+)2 



xi 



a;2 -2/2 



(2/2q+ - 1+)^ + ml 



2/2 



Da = Da 



(18) 



iV^i 



Na, 



Nas 



2/2 -2:2 



2^2 -2/2 



-(Z?2+i?4), 
-(Z?2+i?4), 



2^2 -2/2 



-{D2-Di). 



(20) 



Hard scattering amplitudes for the helicity A -I- A(=0 or 
±1) components for the diagrams Bi{i = 1,2,3) are given 
by 



T, 



(A+A) 

Bi 

(A+A) 



.(A+A), 



(A+A) 



= 7^2 (a; ^ 2/,k+ ^ -1+), 
= TX',^^ (x ^ 2/,k+ ^ -1+), 
= rJA+A)^^ ^ y,k+ ^ -1+). 



(21) 



If one includes the higher twist effects such as intrinsic 
transverse momenta and the quark masses, the LF gauge 
part proportional to 1/fc^ leads to a singularity although 
the Feynman gauge part g^j/ gives the regular amplitude. 
This is due to the gauge-invariant structure of the ampli- 
tudes. The covariant derivative D^ = d^ + igA^ makes 
both the intrinsic transverse momenta, k+ and 1+, and 
the transverse gauge degree of freedom gA± be of the 
same order, indicating the need of the higher Fock state 
contributions to ensure the gauge invariance. However, 
we can show that the sum of six diagrams for the LF 
gauge part(l/fc+ terms) vanishes in the limit that the 
LF energy differences A^, and Aj, go to zero, where A^; 
and Aj, are given by 



A^ = M^ - 



ki 



Xl 



X2 



^ = M2 - M2, 



M-" 



li 



M^ 



2/1 



2/2 



M^. (22) 



In the Appendix, we briefly summarize the proof. 

In this work, we calculate the higher twist effects 
in the limit of A^ = Ay = to avoid the involve- 
ment of the higher Fock state contributions. Our limit 



Aa; — Ay = 0(but ■\/(k^) — /3 ^ 0) may be considered 
as a zeroth order approximation in the expansion of a 
scattering amplitude. That is, the scattering ampHtude 
Th may be expanded in terms of LF energy difference 
A as Tff = [Th](°) + A[ThY^^ + A^[ThY'^^ + • • • , where 
[Th] corresponds to the amphtude in the zeroth order 
of A. This approximation should be distinguished from 
the zero-binding(or peaking) approximation that corre- 
sponds to M = TOi + 7712 and k^ ~ {3 — 0. The point of 



this distinction is to note that [Tff]'-'^-' includes the bind- 
ing energy effect(i.e. kx,l^ 7^ 0) that was neglected in 
the peaking approximation. 

In zeroth order of A^: and Aj,, the net contribution 
from the LF gauge part(l/fc+ terms) vanishes[see Ap- 
pendix] and we only need to compute the Feynman 
gauge(g'"') part, i.e. Na and Nb, for the PQCD analysis 
of meson form factor. The contribution of the Feynman 
gauge to the diagrams Ai is given by 



T 



(A+A) 



J 



y-y(A+A) ^ j^j(X+X) f 6{y2 - X2) 
i-[ ^' ^ \ V2-X2 



N 



(A+A) 



(j/2 - X2)DiD2 



1 
D1D2 



d{x2 - V2) 



•«2 - y2 



0{y2 - X2) + e{x2 - y2){l + 



D12'' 



D.Da D.De 



Axy X _l 



D2 



(1 



D12' 



(23) 



where D12 = Di — D2 and A^y = Ax + Aj,. In the zeroth 
order of A^ and Aj,(i.e. A^; — Ay = 0), Eq. (23) reduces 
to 



T 



(A+A) 



N 



{\+\) 



(2/2 - X2)DiD2 

Similarly, we obtain for the diagrams Bi as 



T, 



(A+A) 



= T 



(A+A) 



{x ^ y, ki 



1| 



(24) 



(25) 



The hard scattering amplitude for each helicity is sum- 
marized as follows: 



^(0) _ 

^H — 



N 



(0) 



(2/2 - X2)DiD2 

8 X2yiy2{x2q\ 



?y, kj 



li 



+ (.T 

+ 2k+ • q+) + X2{xiyi + 2:2^2)1+ • q+ + xiyiml + X2y2'ml + mim2{yi - xi){y2 - X2) 



yi 

+ {x 



^H ^ 



N 



(±) 



(a:2qi 



2k+ • q^)[x2y2{x2(ll_ + 2k+ • q+) - 2x^y2l+ • q+ + (^2 - X2)^Tnl 



+ {x ^^ y, k+ ^^ — Ij 



(y2 -X2)DiD2 

_8 .T2I+ ■ q+ + X2y2ml + xiyiml + mim2{yi - xi){y2 - X2) 

yi (a;2qi + 2k+ • Ci±_)[x2yl{x2(^^ + 2k+ • q+) - 2x^^21+ • q+ + (y2 - X2fml] 



+ {x^y, k_L 



-1±),(26) 



where we neglect the terms such as k^/q^, 1^/q^, and 
kj^ • Ij^/q^ both in the energy denominators and the nu- 
merators due to the fact that k^^ ^ q^ and 1^ <C q^ 
in large momentum transfer region where PQCD is ap- 
plicable [22]. In the hard scattering amplitudes given 
by Eq.(26), the time-ordered 9 function disappears via 
6{x — y) + 0{y — x) = I and there is no singularity in 
timelike region. We also note that the helicity flip con- 
tributions, i.e. mi?7i2(yi — Xi){y2 — 2:2), in the numer- 



ators and the mass terms (^2 — X2)'^m'^ in the denom- 
inators in Eq. (26) give negligible contributions. One 

can easily find that our result Tj^ in the leading twist 
limit reproduces the usual leading twist PQCD result, 
i.e. r^*) = 16/(a;22/2Q'). 



IV. TIMELIKE FORM FACTOR OF A HEAVY 
PSEUSCALAR MESON 



We now consider the PQCD analysis of the timehke 
form factor for the process of e+e~ annihilations into 
two pesudoscalar mesons. The hard contribution to the 
timelike form factor for the electron-positron annihila- 
tions into two pseudoscalar mesons, i.e. e^e" 



MM, 



is obtained as 



Fniq^) = ei/((7^,mi,m2) + 62/(17^, m2, mi), (27) 
with the amplitude I{q^,mi,m2) given by 

/(g^, mi, 7712) == 




where as is the QCD running coupling constant and 
Cf{= 4/3) is the color factor. We note that all the in- 
variant masses in TZq and TZ±i in the spin-orbit wave 
function are replaced by the physical meson mass to be 
self-consistent with the result for the hard scattering am- 
plitude in zeroth order of the LF energy differences A^, 
and Ay, i.e. Ax,Ay — > 0. Again, our analysis should 
be clearly distinguished from the peaking approximation 
(zero-binding limit or /3 ^ 0) which leads to the DA 
(J)m{x) defined by Eq. (9) as the (5- type function with 
Xi — nii/M . In our zeroth order approximation of LF 
energy differences, we consider /3 ^ (i.e. the effect 
of the transverse size of the meson). In the next sec- 
tion (Sec.V), we take the (i values determined from the 
analysis of the mass spectroscopy using our LFQM with 
the variational principle for the QCD-motivated effective 
Hamiltonian [16]. 

Since we neglect the k^ , 1^ and k^ • 1^ terms compared 
to large q^ as we stated below in Eq. (26), we have only 
k±(lj_) • q± and q^ in both numerator and denominator 
terms in Eq. (26). Thus, for convenience in our numer- 
ical calculation, we change the denominator in Eq. (26) 
to include only even powers of kj_ and 1_l- Then the 
numerators of the hard scattering amplitudes T^ and 
T^ in Eq. (26) include both even and odd powers of k^ 
and Ix via the terms (k^ • qx)™ and (1^ • q±)™, where 
m ^integer. Schematically the hard scattering amplitude 
is then given by 



Th^ 



JV((k^-q^)"^,(L-q^)") 
i?(qi",(ki.qi)2«,(U.q^)2n)' 



(29) 



where n and m are integer. In Eq. (29) we show only es- 
sential terms, k^ • qx and k^ • qx in the numerator to ex- 
plain how to obtain the nonvanishing contributions from 
the ordinary and the higher helicity contributions. That 
is, as one can see from Eqs. (13) and (14), by combining 

T^ \T^ ] in Eq. (29) with the spin-orbit wave function 
7^o[7^±] in Eq. (14) to get Th in Eq. (13) or Eq. (28), 
Tfj \Th ] should have even [odd] powers of (k^ • q±) 



and (1^ • q^) in the numerator since 7^o[^±] includes 
even [odd] powers of kj^ and 1^. As a result, we can get 
Th in Eq. (28) as a function of even powers of k_L and 
1^ for both ordinary and higher helicity contributions. 
We then analytically continue to the timelike region by 
changing qx to iqj^(or q"^ -^ —q^ in this case) in the 
form factor. 

In terms of FM{q^) given by Eq. (27), the cross section 
of the pseudoscalar meson pair(MM) production in the 
unpolarized e+e~ annihilations is given by 



da 

dn 



(e+e" -^ MM) 



3^ 
327r 



'^e+e--*fi+fi- Sm' 



,2m2 



\FM{q')\\ 
(30) 



where P 



^/l - 4M2/g2 and CTe+e--»^+^- = 
with Sbcam = y?/2. 

In the peaking approximation, the transverse momenta 
of the quark and antiquark are neglected and the longitu- 
dinal momentum fractions are given by xi = yi = mi/M 
and X2 = y2 ~ m2/M with M = mi + m,2. In this ap- 
proximation, the higher helicity contribution to the hard 
scattering amplitude also vanishes and the ordinary he- 
lictiy contribution to the hard scattering ampitude given 
by Eq. (26) can be rewritten as 



[71(0)1 

[^ H J peaking 



8 ~xlyiy2q^ 



xiyiml +X2y2m\ 



xlyiyl 

+ {x ^ y) 
32M7 / M 
\m2 



g" 



q^ 2m2 



4M2 mi 



(31) 



where 7 — m.im,2/M. 

Therefore, the peaking approximation of the timelike 
form factor of a heavy pseudoscalar meson is given by 

[fM]pcaking(9^) OC ei dxdyS{xi - mi/M)[TP]pc:aking 



X(5(y, -m,/M) + (1^2) 



OC l^ei 



— 1- 



+ 62 



M 

mi 



M'' 


4 
\ 


m2j 


^ 4 






1- 


L/ 





2m2 



4M2 mi 



2?TT,l 



4M2 



m2 



(32) 



This reproduces the result obtained by Brodsky and 
Ji [11]. The form factor zero in this approximation occurs 
at 



nil 
2m2 



62. 

ei J 2mf 



4M2 



1+ ( ^ , . 



(33) 



Even though m^i > ?7i2, the 62 contribution is not negligi- 
ble for the heavy-heavy pseudoscalar meson system such 
as Be- The reason for this is because the timelike form 
factor of heavy pseudoscalar meson encounters a zero and 
62 contribution in the region near the form factor zero 



has non-negligible effect. However, for the heavy-light 
quark system such as B and D mesons, the light quark 
contribution(i.e. 62) can be safely neglected. 



V. NUMERICAL RESULTS 

In our numerical calculations, we use the model pa- 
rameters {fTT-QQ, (3qq) for the linear confining potential 
obtained from the meson spectroscopy with the varia- 
tional principle in our LFQM [16]. Our model parame- 
ters are summarized in Table I. As mentioned earlier, we 
should note that the r.m.s value of the transverse mo- 
mentum in our LFQM is equal to the gaussian (3 value. 



(ki)QQ = 

The shape of the quark DA which depends on the /3 
value is important to the calculation of the cross section 
for the heavy meson pair production in e^e~ annihila- 
tions. We show in Fig. 2 the normalized quark DA of 
B and Ds[D] mesons with different values of (3. For 
the quark DA of _B-meson, we compare our LFQM re- 
sult (solid line) with the small /3 value result close to the 
peaking approximation, e.g. /3 = 0.1 GeV(dashed line). 
As one can see, our LFQM result for the quark DA of B 
meson shows sizable deviations from the peaking approx- 
imation. For the quark DA of D^ [solid line] and D [dotted 
line] mesons, the peak for (?!)_d(x) is located to the right of 
that for 4)o^{x). This indicates that the c-quark carries 
more longitudinal momentum fraction in D than in Dg 
as one may expect. We also show the ratio of (pjj^ (x) and 
(/)£) (x) [dash-dotted line] for the sake of comparison. 

In Fig. 3, we present also the normalized DA of vari- 
ous ^Sq quarkonium((7q) states obtained from our LFQM 
parameters given in Table I. To explain the discrep- 
ancy between the NRQCD prediction [5] and the Belle 
measurement [2] for a{e^e~ — > J /tp + rjc), Bondar and 
Chernyak(BC) [14] reasoned that the discrepancy may 
be due to the extreme (5-function-like charmonium DA 
adopted from NRQCD and claimed that they can fit the 
Belle data by choosing rather broad DA for the char- 
monium state. Interestingly, our LFQM prediction for 
(pec shown in Fig. 3 looks quite similar to BC's result 
in a sense that the DAs for heavy quarkonium states 
differ from the (5-function-likc DA. In our model calcula- 
tion, the DA gets narrower as /3 gets smaller. Also the 
timelikc form factor i^M(<Z^) with small (3 value decreases 
faster than that with large j3 value. Since the cross sec- 
tion (ye+e---yMM is proportioual to |FM(g^)|^, the cross 
section with small /3 is small compare to that with large 
/3. Thus, the cross section for e+e^ — > MM can be in 
principle enhanced by broadening the quark DA. 

In Fig. 4, we compare the results of (j)cc{x) in more 
detail. The solid and dashed(dash-dotted) lines repre- 
sent our LFQM result and BC's result [14] with u^ = 
0.15(u^ = 0.3), respectively, where the parameter ■(;[14] 
represents the characteristic quark velocity in the bound 
state. Although there is a similarity in the quark DA be- 
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FIG. 2: Normalized quark distribution amplitudes of _B[top] 
and Ds [D] [bottom] mesons for different values of the gaussian 
parameter /3. 



tween ours and BC's results(in particular, their u^ — 0.15 
result), there is rather substantial difference near the end- 
point region between ours and BC's results. Since the 
PQCD hard scattering amplitude is typically very sen- 
sitive to the end-point values of DA, it may not be so 
difficult to imagine that BC's prediction of double charm 
production cross sections would have made a difference 



TABLE I: The constituent quark masses mg[GeV] and the gaussian paramters l3{- 
obtained from the variational principle, q—u and d. 



\/{k^)[GeV] for the linear potential 
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FIG. 3: Normalized quark distribution amplitudes of qq with 
our LFQM parameters Pqq. 



FIG. 4; Normalized quark distribution amplitude of cc(solid 
line) compare to those(dashed, dash-dotted lines) obtained 
from Bondar and Chernyak [14]. 



depending on what w^ value they have used. To fit the 
Belle measurement [2] for (T(e+e^ -^ J/tp+rjc), they used 
v"^ = 0.3 rather than v"^ — 0.15. However, as we stressed 
earlier, the way that BC[14] handled the transverse mo- 
mentum effect is different from ours since they integrated 
out the transverse momentum seperately in the wave 
function part and in the hard scattering part while we 
didn't factorize the hard and soft parts but integrated out 
the transverse momentum for the whole amplitude. We 
think that a consistent analysis with (3^0 (or non-zero 
binding energy) should follow our non-factorized formu- 
lation (see e.g. Eq.(TF)). 

In Fig. 5, we show the cross sections for the exclu- 
sive Df D^ and D^D^ pair productions in e+e^ anni- 
hilations. The solid line represents the results including 
both ordinary and higher helicity contributions while the 
dotted line corresponds to the result of the ordinary he- 
licity contribution only. The dashed line for a^e'^e^ -^ 
D~^D~) represents the lower limit(i.e. f3 — 0.16 GeV) 
for the form factor zero to occur above -^s = 10.6 
GeV. The small black circle for cr(e+e^ -^ D^D^) 
represents the upper limit obtained from Belle [3], i.e. 
crExp.(e+e~ -^ D+D-) < 0.04. The higher helicity 
contribution to the cross section, i.e. the difference 
between solid and dotted line, is more pronounced in 



e+e^ -^ D^D^ than in e+e^ -^ D^Dj^ especially 
near the turnover point (or the form factor zero point). 
In general, the higher helicity contribution to the me- 
son form factor increases(decreases) as quark mass de- 
creases (increases). For instance, while the higher helicity 
contribution to the hard scattering amplitude is negligi- 
ble in the B meson form factor, it is not negilible in the 
pion form factor. However, the most significant in our 
analysis is the transverse momentum effect which delays 
the turnover point (see Eq. (33) for the peaking approx- 
imation). For instance, the turnover for Ds[D] meson 
occurs near 6. 3 [11. 3] by going beyond the peaking ap- 
proximation while the corresponding turnover point is 
near y/s/M ^ 2.8[4.03] for the peaking approximation. 

Numerically, we obtain the cross sections for D^Dj 
and D~^D~ pair productions at ^/s = 10.6 GeV with our 
LFQM paramters(i.e. /? = 0.5016 GeV for Ds and 0.4679 
GeV for D) as follows 



<j{e+e- ^ DtD-) = (S.Ot^;^) x 10-4[pb] 
a{e+e- -f D+D-) = (0.02 ± 0.01)[pb], 



(34) 



for the strong coupling constant as = 0.2 ± 0.05. Simi- 
lar Us values were used in Refs. [23, 24]. Our result for 
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FIG. 6: Lower bound for the shape of ijf>_D (a;) [dashed hne] de- 
duced from the assumption of a{e'^e~ -^ D'^ D~) < 0.04[pb] 
and (da/d^/s) < at ^s = 10.6. 



__(x,k.j^): p=0.4679 GeV " 
(orcKnary+higher helieities) 
rLFnM(''.V-P=0-4679GeV 
{ominary helicity only) 

,Q„(x,l4):P = 0.16 GeV 
(dftnnary+higher helieities) 
Belle Exp. < 0.04(upper limit) 




FIG. 5: Cross sections for DfDg [top] and D^D [bottom] 
pair productions in e+e^ annihilations with Os — 0.2. 



CT(e+e^ -^ D^D^) is consistent with the recent experi- 
mental data from Belle, crExp.(e~''e~ -^ D^D^) < 0.04. 
Since a{e~^e~ -^ D^D^) gets larger as beta grows, the 
upper bound of (TExp. from Belle provides a constraint on 
the maximum /3 value modulo the dependence on ag. 

If the cross section for the D meson pair produc- 
tion satisfies cr(e+e^ -^ D^D^) < 0.04 and its slope 
with respect to the momentum transfer is negative, i.e. 



(da/dy/s) < at y/s = 10.6 GeV, then we could also set 
the lower bound for /3 value as /3 > 0.16. The shape of 
the D meson quark DA corresponding to (3 — 0.16 GeV 
is shown by dashed line in Fig. 6. If /3 < 0.16 GeV, the D 
meson quark DA approaches to the (5-type function but 
(da/dy/s) > at y/s — 10.6 GeV due to the form factor 
zero occuring at y/s < 10.6 GeV. Due to the occurence of 
form factor zero for the heavy pseudoscalar meson pair 
production [4, 11, 25], more experimental data around 
y^ — 10.6 GeV are necessary to check the slope of the 
cross section. More data around ^/s = 10.6 GeV would 
further constrain the shape of D meson quark DA. 

How about B mesons? We should point out that the 
PQCD result of the cross section for B^B~ pair produc- 
tion may not be trustworthy because y/s — 10.6 GeV is 
too close to the threshold energy of B^B^ pair produc- 
tion. As expected, the gluon momentum transfer from 
the heavy quark to the light quark in B^ B^ pair produc- 
tion at ^/s — 10.6 GeV turns out to be only around a few 
hundred MeV close to the scale of Aqcd- On the other 
hand, the gluon momentum transfer in D7,,D7, pair pro- 
duction at y/s = 10.6 GeV is much larger than the scale 
of Aqcd- By going beyond the peaking approximation, 
the average gluon momentum transfer gets even larger 
due to the transverse momentum effect. This may jus- 
tify our PQCD analysis for D'^'^.DT', pair production at 

Vi = 10.6 GeV. 

Although the absolute value of the cross section for 
i3[i3s]-meson may not be reliable near y/s — 10.6 GeV, 
it seems interesting to discuss the behavior of the ratios 
of cross sections such as (T(e+e^ -^ B^ B'^) / a{e^ e'^ -^ 
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FIG. 7; Our predictions [closed symbols] on the cross section 
ratios for various heavy pseudoscalar ineson(_B, Bs, i3,and Ds) 
pair productions compared to the peaking approximation re- 
sults[open symbols] near y'i = 10.6 GeV. 
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B+B- 



B+B-) and a{e+e- 
In Fig. 7, wc show our predictions [closed symbols] on 
the cross section ratios for various heavy pseudoscalar 
meson ( i?, i3s, Z), and Ds) pair productions in e+e~ an- 
nihilations near ^/s = 10.6 GcV, i.e.y^ ranging from 
10 to 12 GeV and compare our results with those[open 
symbols] obtained from the peaking approximation [11]. 
The open and closed diamond symbols for a{e'^e^ -^ 
B'^B'^)/a{e^e^ — > B^B~) are on top of each other and 
their values are almost equal to 1. In fact, the cross sec- 
tion ratios involving light quarks u and d such as {D^ , D) 
and (-6°, B) cases are close to 1. This is due to the negli- 
gible contribution from the diagrams where the photon is 
attached to the light quarks. As the replacement of light 
quarks by the strange quark makes those diagrams non- 
negligible, the cross section ratios for the cases of {Ds^ D) 
and (Bs,B) deviate from 1 appreciably. However, most 
significant is again the transverse momentum effect which 
is pronounced in the case of D[Ds] meson pair produc- 
tions compare to B[Bs] meson pair productions. In par- 
ticular, the deviation between the open and closed sym- 
bols for the case {Ds,D) is quite dramatic compare to 
the case of {Bs,B). 



found to be relat ed to the transverse momentum via 
Pqq = w (k^)QQ. This relation naturally explains the 

zero-binding energy limit for the zero transverse momen- 
tum, i.e. (Mq ) = (mi -I- 7712)^ and Xi = nii/M for (3 = 0. 

However, the heavy quark DA is sensitive to the value 
of l3 and indeed substantially broad and quite different 
from the J-type DA according to our LFQM based on 
the variational principle for the QCD-motivatcd Hamil- 
tonian [16, 17]. If the quark DA is not an exact 5 func- 
tion, i.e. k^ in the soft bound state LP wave function 
can play a significant role, the factorization theorem is 
no longer applicable. To go beyond the peaking approx- 
imation, the invariant amplitude should be expressed in 
terms of the LF wave function ^(x^, k^^) rather than the 
quark DA. 

In going beyond the peaking approximation, we 
stressed a consistency by keeping the transverse momen- 
tum k^ both in the wavefunction part and the hard scat- 
tering part together before doing any integration in the 
amplitude. Such nonfactorized analysis should be distin- 
guished from the factorized analysis where the transverse 
momenta are seperately integrated out in the wavefunc- 
tion part and in the hard scattering part. Even if the used 
LF wavefunctions lead to the similar shapes of DAs, pre- 
dictions for the cross sections of heavy meson productions 
would apparently be different between the factorized and 
nonfactorized analyses. 

In this work, we compared our non-factorized analysis 
with the usual factorized analysis based on the peaking 
approximation and found substantial difference between 
the two in the calculation of the cross section for the 
heavy meson pair production. We also discussed the 
higher hclicity contribution to the cross section. Our 
analysis provided a constraint on the size of quark trans- 
verse momentum inside the D meson from the recent 
Belle data, CTExp.(e+e" -^ D+D') < 0.04. More ex- 
perimental data around ^/s = 10.6 GeV would further 
constrain the shape of D meson quark DA and test our 
LFQM prediction. 
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VI. SUMMARY AND CONCLUSION 

We investigated the transverse momentum effect on 
the exclusive heavy meson pair productions in e~^e~ 
annihilations within the framework of LF PQCD. The 
gaussian parameter (3 in our model wave function is 



APPENDIX A: PROOF OF VANISHING 
LIGHT-FRONT GAUGE PART IN M = Mo LIMIT 

The contribution of the LF gauge part(l/A:+ term) to 

the hard scattering amplitude for diagrams A = X]i=i ^i 
obtained from Eqs. (17) and (20) is given by 
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In terms of the LF energy differences, the relevant energy 
denominators can be rewritten as 
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Eq. (Al) and the corresponding LF gauge part to the 
diagrams Bi lead to singularties. Fortunately, however, 
in the zeroth order of A^, and Aj,, i.e. Aj; = Aj, -^ 
limit, one can see that the energy denominator term in 
0{x2~y2) in Eq. (Al) vanishes, which leads to 9{y2^X2) + 
9{x2 — 2/2) = 1- Thus, the LF gauge part contribution to 
the diagrams A — X]i=i ^i becomes 



T 



(i/fc+) 



1 



(2/2 - X2Y D2 



(A3) 



and similarly we obtain 



n 



for the diagrams B 
tion 1)2 + -Dg = Aj, ^ 



(2/2 - X2y Dq' 



(A4) 



= X]i=i ^i ■ Finally, from the rcla- 
Aj, , one can see that the net contri- 

(i/fe+) , (i/fc+) 



bution from the LF gauge parts, i.e., T^ 
vanishes exactly in the limit of A^, ~ Ay ~ 
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